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Abstract 

Extensive work has been done on the group classification of systems of equations 
in the literature. This paper identifies the gap in the literature which concerns the 
group classification of systems of two nonlinear second-order ordinary differential 
equations. We provide a complete group classification of systems of two ordinary 
differential equations of the form, y” = F(y), which occur in many physical appli¬ 
cations using two approaches which form the essence of this paper. 
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1 Introduction 


Systems of autonomous nonlinear second-order ordinary differential equa¬ 
tions, where the independent variable, usually assumed to be time, does 
not appear on the right hand side of the system, arise in various physical 
problems. This effectively assumes that the laws of nature which hold true 
in the present are presumably applicable in the past and the future. Hence, 
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like all other systems of differential equations, the study of their symme¬ 
try structure poses an important role as their presence in a system allows 
one to reduce the order of the studied equations and also to hud general 
solutions in quadratures. 


Group classihcation studies, dating more than a century back, were hrst 
initiated by the founder of symmetry analysis, Sophus Lie pll2ll3ll^ . These 
studies were long forgotten until Ovsiannikov [ jSIIB] revived the work hve 
decades ago. Lie’s works put emphasis on tackling the group classihcation 
in two ways: 1) the direct way and 2) the indirect way also known as the 
algebraic approach. The direct way involves directly hnding solutions of 
the determining equations and allows one to study all possible admitted 
Lie algebras without omission. On the other hand the indirect way in¬ 
volves solving the determining equations up to hnding relations between 
constants dehning admitted generators. The algebraic approach, as in the 
studies |71l8ll9llT0]G I . takes into account the algebraic properties of an admit¬ 
ted Lie group and the knowledge of the algebraic structure of the admitted 
Lie algebras in order to allow signihcant simplihcation of the group classih¬ 
cation. In one of Lie’s works [I], he gave a complete group classihcation of a 
single second-order ordinary diherential equation of the form y" = /(x, y). 
Later on Ovsiannikov IH] did this group classihcation in a diherent way. 
The method he used, now also known as the direct approach, involved a 
two-step technique where the determining equations were hrst simplihed 
through exploiting equivalence transformations and later on solved for the 
reduced cases of the generators. The same technique was used in a study 
conducted in |tI2] to classify a more general case of equations of the form 
y" = P 3 (x, y] y')^ where P‘^{x^ y; y') is a polynomial of a third degree with 
respect to the hrst-order derivative y'. Observe that sometimes difficulties 
arise in using the direct approach. Sometimes it is difficult to select or 
tease out equivalent cases with respect to equivalence transformations. As 
is observed in the classihcation of a general scalar second-order ordinary 
diherential equation of the form y" = f{x,y;y'), the application of the 
direct technique gives rise to overwhelming difficulties. In this study both 
the direct and indirect techniques are employed. 


Apart from dealing with classihcation problems there is a signihcant amount 
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of research that deals with the dimension and structure of symmetry alge¬ 
bras of linearizable ordinary differential equations [ 171113014015016017] . This 
is also of importance since some nonlinear equations appear in disguised 
forms. 


In addition to extensive studies on properties of scalar second-order or¬ 
dinary differential equations, there are also several researchers committed 
to studying systems of two linear second-order ordinary differential equa¬ 
tions P!H019017020021022023024025| . Surprisingly, the group classihcation of 
systems of two nonlinear second-order ordinary differential equations has 
not yet been exhausted, in particular, the group classihcation of systems 
of two autonomous nonlinear second-order ordinary differential equations 
is not yet complete. Hence this paper considers the group classihcation 
of systems of two autonomous nonlinear second-order ordinary diherential 
equations of the form 


y" = F(y). 


( 1 ) 


The system studied here is a generalization of Lie’s study [2]. Studied cases 
such as systems of two linear second-order ordinary diherential equations 
and the degenerate case which is equivalent to the following 

y” = F{x,y,z), z' = 0 (2) 

are omitted from this paper. We call systems that are equivalent to these 
cases as reducible systems and irreducible otherwise. 

The paper is organized as follows. A preliminary study of systems of two 
nonlinear second-order ordinary diherential equations is tackled hrst and 
is followed by the subsequent group classihcation applied to autonomous 
systems ([I]) of two second-order ordinary diherential equations. The group 
classihcation is divided into two parts depending on the coefficient of the 
inhnitesimal generator. The direct approach is applied on one case while a 
combination of the optimal system of subalgebras and direct approach is 
applied to the other case. The latter part of the paper lists the diherent 
cases with their respective results and is then followed by the conclusion. 
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2 Background study of systems of the form y" = F(a;,y) 


This section focuses on systems of two nonlinear second-order ordinary 
differential equations of the form 


y" = F(^,y) 


(3) 


where 




y = 



^G{x,y,z) ^ 


2.1 Equivalence transformations 

System has the following equivalence transformations: 

(1) a linear change of the dependent variables y = Py with constant non¬ 
singular 2x2 matrix P; 

(2) the change y = y + (j){x) and 5 = 2 ; + and 

(3) the transformation related with the change x = 0(x), y = yi{j{x), 
z = zif{x)., where the functions 0(a;) and ifi{x) satisfy the condition 



2.2 Determining equations 


The determining equations in matrix form for irreducible systems of the 
form (E]) are given by 


2CF, + 3C'F + (((>1 + ('E)y + C) ■ V)F - AY = {"'y + C", (4) 


where the matrix A = (ajj) is constant. The associated inhnitesimal gen¬ 


erator is 


X = 2?(x)cl, + (Ay + C(x))-V, 


where V = {dy, d^Y and ” • ” means the scalar product b • V = bidy^, where 
the summation with respect to the repeated index is used [[22] • 
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The equivalence transformation ([T]) with linear change y = Py, when ap¬ 
plied to equations (ED, reduces equations (HD and its associated inhnitesimal 
generator to the same form with the matrix A and the vector ( changed. 
Equations (ED become 

y" = F(x,y) 


with 


F(x,y) = PF(x,Fy), 

and the partial derivatives with respect to the variables y are also changed 
as follows 


b • V = (Pb) • V. 


Consequently, the determining equations (ED become 


2CF, + a^'F + (((i + eE)y + C) ■ V)F - iF - ry - C" = o, 


where 


A = PAP-\ C = PC 

and the associated inhnitesimal generator is also changed as follows 


X = 2C(;r)a,+ (iy + C(x))-V. 


As in [ E2IE5] and in the succeeding pages, this transformation places a very 
important role in the group classihcation process. 

From the study 122], the systems of two nonlinear second-order ordinary 
differential equations are equivalent to one of the following 10 types listed 
below. Looking closely at these systems, there is a necessity to conduct an 
initial study where the systems of two equations do not depend on x. This 
forms the core of this paper. 
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F and G 


Relations and conditions 


u = ye V = ze 


1. F = e“^/(«, n), 

G = e^^g{u, v) 

2. F = e^^{cos{cx)f{u,v) + sm{cx)g{u,v)), 

G = e“*(— sin(cx)/(u, v) + cos{cx)g{u, v)) v = e~^^{y sin(cx) + z cos(cx)), 

a, c 7 ^ 0 are constant 


a, b are constant 

(ycos(cx) — 2 ;sin(cx)), 


3. F = e'^^{f{u,v) + xg{u,v)), 

G = e°'^g{u, v) 

4. F = {y + hi{x))f{x,v) - h'l, 

G = (z + h 2 {x))g{x, v) - /i 2 

5. F = {y + hi{x))f{x,v) - Fl{x), 

G = h 2 {x) In (y + hi{x)) + g{x,v) 


6 . 


F = 


G = 


Kjx) 

hi{x) 

K{x) 

hi{x) 


y + fix,v), 
y + g{x,v) 


7. F = e““(cos(cii)/(x, v) + sin(c«)y(x, v)) 
{cu)f{x, v) + cos 
-f{x,v)+g{x,v) 


G = e““(— sin(cii)/(x, v) + cos{cu)g{x, v)) a, c / 0 are constant 

y 


F = - 

z + hi{x)' 

G = -h'lix) + f{x,v) 


9. 


F = 


K{x) 2 


+ uf{x,v) +g{x,v), 


u = e °'^{y — zx), V = ze 
a is constant 

V = {z + h2{x)){y + /ii(x))“, 
a 7^ 0 is constant 

V = z — h2{x) In {y + hi{x)) 

V = z- hi{x) / 0 

hi[x) 

y = sin(cii), z = cos(cu), 
a, c 7^ 0 are cc 

V = z + hi (x) 

(z + hi{x)Y 


z + hiix) 

u= , . , v = y 

h2[x) 

/i2(x) 7^ 0 


h2{x) 


y = uve , z = ve 


G = —h"{x) + h 2 {x)u + /(x, v) 

10. F = e^{uf{x, v) + y(x, x)), 

G = e“/(x, v) 

For all the cases, hi, / 12 , / and g are arbitrary functions of their arguments. 
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3 Autonomous systems ([T]) of two nonlinear second-order ordinary differen¬ 
tial equations and their group classification 


Since for autonomous systems = 0, the determining equations for au¬ 
tonomous systems have the form 

Se'F + ({(/! + i'E)y + C) ■ V)F - AF - f'y - C" = 0. (5) 

This also implies that the generator is admitted by system ([T]). 

Differentiating the determining equations (E]) with respect to x, the group 
classification study is reduced into two cases, namely, 

( 1 ) the case with at least one admitted generator with ^ 0] and 

( 2 ) the case where all admitted generators have = 0 . 

For the hrst case, the direct approach by Lie is utilized, whereas for the 
second case, a combination of the optimal system of subalgebras of the Lie 
algebra and the direct method is used. 


3.1 Systems admitting at least one admitted generator with 7 ^ 0 


For the case with at least one generator, 7 ^ 0, we initially consider the 
differentiated determining equations (Ej) with respect to x and divide them 
by The determining equations become 

»F+((,4).v)r-f,-C.o, (.) 

Fixing X, and shifting y and z, equations (Ej) are reduced to 

3F + (y • V)F — ay — b = 0 
where vector b = ( 6 , c)^, with a, 6 , c constant. 
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The general solution of these equations is 


F = l + J + 

G = I + ^ + z-^giu), 

where u = z/y and f'g' 7 ^ 0. It is easy to see that if f'g' = 0, the studied 
system is equivalent to a reducible case. The functions F and G are then 
substituted into the determining equations (E]). The determining equations 
are then solved directly in order to hnd generators admitted by equations 
([1]). The first part of the determining equations is given as follows: 

i'" - = 0 , ( 8 ) 



{Ciu - C 2 ) fu + 3 Ci/ = 0 , 

{u^Ci - uC, 2 )gu + K 2 g = 0 , 

VlC,i — 126^' — 3a(i + 4aii6 + 4ai2C = 0, 
I 2 C 2 - 12c^' - 3aC2 + 4a2i6 + 4a22C = 0, 


(9) 

( 10 ) 


{aiiu'^ + ai2M^ — a2iw^ 
[aiiu^ + ai 2 W^ — a2iu 


a22U^)fu + (4aiiM^ + ‘iai2U^)f + ai2^ = 0, 

a22U^)gu + a2\u^f + (3021 + ^a22u)g = 0 . 


( 11 ) 


From equation (jS]), it can be seen that the general solution of ^ depends on 
three values of a, i.e., a = 0, a = —p^ and a = where p 7 ^ 0. For a = 0, 
the general solution of ^ is 


^ — ^2x‘^ + 

where ,^2 7 ^ 0, ,^ 1 , ^0 are constant. For a = —p^, the general solution of ^ is 

^ = 6 cos(px) + <^2 sin(pa;) + <^o, 

where ^2 7 ^ 0, 7 ^ 0, ^0 are constant. Lastly, for a = p^, the general 

solution of ^ is 

^ + 6 , 
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where ^2 7 ^ 0, 7 ^ 0, ^0 are constant. Subsequently the determining equa¬ 

tions (Ej) lead to the study of two cases where: ( 1 ) there exists a generator 
with Cl 7 ^ 0 and ( 2 ) all generators have Ci = 0- 

Considering the case where there exists a generator with Ci 7 ^ 0, we di¬ 
vide by Cl and differentiate the equations (El) with respect to x to obtain, 
C2 = ^Ci? where /c is a constant. Substituting this back to equations (E]), 
one obtains / = /o(w — k)~^ and g = gou^{u — k)~^. Also, differentiating 
equations (HUI) with respect to a;, it follows that c = kb. From here, one can 
verify that this is a reducible case. 


Consider that all generators have Ci = 0 . From equations (E]), it follows 
that C 2 = 0. Differentiating equations (pUj) with respect to a;, it immediately 
follows that 6 = c = 0. From equations (1111) . the equivalence transformation 
y = Py^ where P is a constant nonsingular 2 x 2 matrix, is utilized in order 
to obtain the general solution of / and g. Note that the constant matrix 
A is reduced to one of the following real-valued Jordan forms 


Ji 


f n ^ 


( 1 ^ 


( 1 \ 

ail 0 


ail 1 


ail 1 


5 J 2 = 


5 Ji = 


CM 

CM 

0 


^-1 an ^ 


^ 0 ail ^ 


( 12 ) 


The general solutions for / and g are listed as follows: 


Jordan form 

f 

9 

Ji 


r _ 

Jiu '^-1 

J 2 

ifoy - fiz)r{y, z) 

{foz + fiy)r{y, z) 

Ji 

e~^{fou-^ + fiu-^) 



In this instance, T{y, z) = y {y‘^ + z‘^) and /o 7^ 0 , /i 7^ 0 , 7 7^ 1, 

(a 7 ^ 1 are constant. 

Excluding reducible systems, the classes of functions F and G of equations 
(Ej) admitting a corresponding Lie group and the extension of the kernel of 
the admitted Lie algebras are obtained as seen in Table H The kernel of 
the admitted Lie algebras consists of the generator Xi = d^, which will be 
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omitted on the list. The extension of the kernel is listed as follows: 

Y 2 = 2 xdj; + ydy + zdz, Y 3 = x{xdx + ydy + zd^), 

Ya = -fydy + zd^, Y 5 = {ay - z)dy Y {y Y az)d^, 

^6 = ( 2 / + 2 xdx — sin 2 x{ydy + zd^), 

Ys = sm2xdx + cos2x{ydy + zd^), Yg = e~‘^^{dx - {ydy + zd^)), 

Yio = e^^{dx + ydy + zd^). 

The Lie algebras >25 ^ 3 , ^ 7 , ^ 8 ? ^ 9 ? and Yig are associated with the coeffi¬ 
cient ^ and the Lie algebras I4, Ys, and Yq are related to the type of Jordan 
form of matrix A. 


3.2 Systems where all admitted generators have = 0 


Note that the action of equivalence transformations coincides with the ac¬ 
tion of group automorphisms. For the direct approach, sometimes it is 
difficult to select out equivalent cases with respect to equivalence transfor¬ 
mations. Fortunately, if the algebraic structure of the admitted Lie algebra 
is known, then using the algebraic approach aids in simplifying the group 
classihcation problem. Thus, for hnding the group classihcation of systems 
of two autonomous nonlinear second-order ordinary differential equations 
with all admitted generators satisfying = 0, the one-dimensional opti¬ 
mal system of one parameter subgroups is utilized and is proceeded by the 
direct approach. Firstly, the two-step algorithm of Ovsiannikov |26l 


IS em¬ 


ployed here, for which the optimal systems of subgroup and group invariant 
solutions are reconstructed. Invariant solutions are then substituted back 
into the determining equations (ED where the direct method is used to hud 
all possible admitted Lie algebras. 


Firstly returning to the analysis of the determining equations (ED, since 
= 0, it follows that ^ = /ci + k 2 X, where ki and kp. ar e constant. This 
property of the coefficient forces ^ to become constantThe determining 


**See in Appendix. 
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equations (E]) are then reduced to 


3*2F + (((^ + A2£)y + k)-V)F-^F = 0 (13) 

with the following admitted generator 

X = 2{ki + k2x)dx + (^y + k) • V (14) 

where the matrix ^4 is a vector k, ki and /c 2 are constant. By rewriting 
(fT4l). the generator can be represented as 

X = E qX, (15) 

i=\ 

where Ci{i = 1,.. .8) are constant. Corresponding to the constants Ci{i = 
1, ...,8), the basis operators of the Lie algebra are as follows: 


Xi = 4 X2 = xd, Xa = dy X4 = d, 
X 5 y^y Xs zd;^ X 7 zOy Xg ydz- 


(16) 


From here, the one-dimensional optimal system of one parameter subgroups 
of the main group of system ([I]) with = 0 is constructed. The commu¬ 
tators of the basis operators are 


[Xl, X2] 

= Xi, 

[X5, X7] = 

[X3, X5] 

= X 3 , 

[X5, Xs] = 

[X3, Xs] 

= X 4 , 

[Xs, X7] = 

[X4, Xs] 

= X 4 , 

[Xs, Xs] = 

[X4, X7] 

= X 3 , 

[X7, Xs] = 


-X7, 

Xg, 

X7, ( 17 ) 

-Xg, 

6 - X 5 . 


The following inner automorphisms 74j(i = 1,..., 8) of the above Lie algebra 
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are found without difficulties: 

Ai : Cl = Cl — aiC2, 

^2 : Cl = e“ 2 ci, 

^3 : C3 = C3 — a3C5, C4 = C4 — a3C8, 

A4 : C3 = C3 — fl4C7, C4 = C4 — U4C6, 

^5 : C 3 = e“5c3, C 7 = e“5c7 cg = e'^^cg, 

lo j 

^40 : C 4 = e®®C 4 , C 7 = e““®C 7 cg = e^^cg, 

^ 4 ? : C3 = C3 + a7C4, C5 = C5 + a7Cg, ce = ce — ajcg, 

07 = 07 — OjOs + a 70 Q — a7C5, 
v 48 : C4 = C4 + agC3, C5 = C5 — agC7, ce = ce + a8C7, 
cs = Cg — a8C7 — agce + ugcs. 

Note that (z = 1, 8 ) are the parameters on which the transformations 
of the group depend on. Apart from these automorphisms, we have the 
following involutions: 

£^1 . z z\ C4 C4, O7 O7 , Cg Cg, 

E 2 -. y = -y\ C 3 = -C 3 , 07 = - 07 , cg = -cg; 

: X = —x\ Cl = —ci; 

E 4 . . y z^ z ^1 C3 C4, C4 C3, C5 C0, C5 C5, 07 Cg, Cg 07 . 

We study the way in which the coefficients of equation (ITN) are changed 
under the action of inner automorphisms of the group above. Here and 
further on, only changeable coordinates of the generator are presented. 
Looking closely at the commutators, the Lie algebra Lg, which is composed 
of the generators Xi {i = 1 ,..., 8 ), can be split into 2 subalgebras L 2 0£6 = 
{Xi, X 2 } 0 {X 3 , X 4 , X 5 , X 0 , X 7 , Xg}. Note also that Lq can be decomposed 
further to £ 40/2 = {X 5 , X 0 , X 7 , Xg} 0 {X 3 ,X 4 }, where £4 makes up a 
4-dimensional subalgebra and I 2 is ideal. 

Let us first study the 4-dimensional subalgebra £4 = {X 5 , X 0 , X 7 , Xg}. We 
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consider this study here due to a misprint found in the classihcation of 
this Lie algebra in p7]. Now consider the operator X of a one parameter 
subgroup of the form 


X — C 5 X 5 + cqXq + C 7 X 7 + cgXg. (19) 

Automorphisms A 5 up to Ag are made use of in order to hud the one¬ 
dimensional optimal system of subalgebras of this Lie algebra. From the 
automorphisms A 5 and As, one can hud the invariant C 7 Cg = C 7 Cg, which 
leads one to consider the following cases: 

(а) cjcg > 0 

(б) C7Cg < 0 

(c) C 7 Cg = 0 . 

Then utilizing the invariant of A 7 and Ag, which is C 5 + ce = C 5 + cs, one 
can obtain relations of C 5 and cq. Upon further computations using auto¬ 
morphisms, it can be verihed that for case (a), the coefficients of equation 
(113) satisfy C 5 — ce 7 ^ 0, C 7 = 0 and cg = 0. For case ( 6 ), it follows that 
C 5 = ce, C 7 = — 1 and cg = 1. For case (c) if C 5 7 ^ cq then C 7 = 0 and 
Cg = 0. If C 5 = ce, then C 7 = 1 and cg = 0. The involutions are also utilized. 
Hence, the following one-dimensional optimal system of subalgebras of the 
Lie algebra L 4 is obtained: 


1 . X 5 + aXQ where — 1 < a < 1 

2. q:(X 5 + Xq) + Xg — X 7 where a > 0 

( 20 ) 

3. /^(Xs + Xe) + X 7 where /I = 0,1 

4. 0. 

Note that the 0 element is considered on this list [22]. There is a necessity 
to include this element on the list as when the direct sum L 40/2 is applied, 
more subalgebras of the Lie algebra Lq may appear on the list. 

Remark: As the action of the above automorphisms coincides with the 
action of the equivalence transformations, it is possible to get the opti¬ 
mal system of one-dimensional subalgebras of the Lie algebra L 4 using 
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the latter. From the determining equations (USD and the utilization of the 
equivalence transformation y = Py, where P is a nonsingular 2 x 2 matrix 
with constant entries, the matrix of coefficients of ( IT^ 


( \ 

C5 C7 

^ Cg Cq ^ 

is reduced to one of its real-valued Jordan forms ( IT^ . Looking closely at 
(120), subalgebra 1. coincides with Jordan matrix ffi, subalgebra 2. coincides 
with Jordan matrix J 2 , and subalgebra 3. coincides with Jordan matrix J 3 . 


3.2.1 Optimal system of subalgebras of the algebra Lq = {X 3 , X 4 , X 5 , Xg, X 7 , Xg} 

After obtaining the one-dimensional optimal system (I2UI) of subalgebras of 
the Lie algebra L 4 = {X 5 , Xs, X 7 , Xg}, the next step is to combine L 4 with 
the ideal h = {X 3 ,X 4 }. Here, again Ovsiannikov’s two-step method | 2 B] 
is applied. Hence, for the study of the one-dimensional subalgebras of the 
Lie algebra Le, the study is reduced to analyzing the following elements: 

1. C3X3 + C4X4 + X5 + aX^ where — 1 < a < 1 

2 . C3X3 + C4X4 + a{X5 + Xq) + Xg - X7 where a > 0 

3. C3X3 + C4X4 + /J(X 5 + Xq) + X7 where /J = 0,1 

4 . C3X3 + C4X4. 


Using automorphisms A 3 and A 4 and the involutions, the list of one¬ 
dimensional subalgebras of the Lie algebra Lq = {X 3 , X 4 , X 5 , Xe, X 7 , Xg} 
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is obtained as follows: 


1. X5 + aXQ where — 1 < a < 1 

2. X4 + X5 

3. Xg - X 7 

4 . /3X3 + a{X5 + Xq) Xg — Xj where (3 = — 1 , 0 , 1 , a > 0 

5 . 13X4^ + X7 where /5 = 0,1 

6. X5 + Xs + X7 

7 . X 3 

8 . 0 . 


Again, it is necessary to study the element 0 of the subalgebras of the Lie 
algebra Lq as this may generate additional elements when Lq is combined 
with L 2 . 


3 . 2.2 Optimal system of subalgebras of the algebra = {Xi, X2, X3, X4, X5, Xg, X7, Xg} 

Combining Lq with L2 and keeping in mind that for autonomous systems 
Xi is already admitted, the following elements comprise the list of one¬ 
dimensional subalgebras of the Lie algebra Lg: 

1 . 7 X 2 + X 5 + aXQ where — 1 < a < 1 

2. 7 X 2 + X 4 + X 5 

3 . 7X2 + Xg — X7 

4. 7X2 + l3Xg + q:(X 5 + Xq) + Xg — X 7 where (3 = —1,0,1, a > 0 

5. 7X2 + ;dX4 + X 7 where /5 = 0,1 

6 . 7X2 + X 5 + Xq + X 7 

7 . 7X2 + X3 

8. X2. 


15 



Using this list of subalgebras, the next step is to obtain invariant solutions 
F and G of the determining equations (HSj). These functions are substituted 
into the determining equations (Ej), which are solved completely in order 
to hnd all other generators admitting equations ([I]). 

3.2.3 Representations of systems of two nonlinear second-order ordinary differential 
equations with all generators having f" = 0 

From ( 1 T 51 ). Ci{i = 1 ,..., 8) are the coefficients of the generator chosen from 
the above list of subalgebras. Only one subalgebra is presented in this paper 
as computations for the other subalgebras are done in a similar way. 


3.2.3.1 Subalgebra 1. with the generator 7 X 2 + X 5 + aX^ where —1 < a < 1. 
For this case, the determining equations ( fT^ become 

yFy + azFz - (27 - 1)F = 0 
yGy + azGz — (27 — a)G = 0 . 

The general solution of these equations is 

F{y, z) = f{u)y^-F and G{y, z) = g{u)y^-F^ 

where u = y^/z. Notice that /' 7^ 0 else it is equivalent to a reducible 
case. Substituting these functions to the determining equations (E]), the 
following initial determining equations are obtained 

y^^auiauf + (1 - 27)/ - ug) 

+y^+^u{{aan + {a - l)^i - 022)^/' - 2(7011 + (7 - 2)^i)/) 
ffy^CMauf + (1 - 27)/) - y^a2iu^f' - yC,2U^f = 0, 
y‘^^ai2{aug' F {a - 2j)g) 

Fy^^^u{{aau F {a - l)^i - a22)ug' - ((a - 27)011 + (a - 27 + 3 )^i - 022)5') 
Fy^uCi{aug' F {a - 2 -f)g) - y‘^a2iu{u^g' F f) - yC,2VZ'g' = 0 
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In order to split these determining equations, one needs to study relations 
between the powers of y. Thus, one needs to evaluate the following cases: 
( 1 ) a = 0 , ( 2 ) a = ( 3 ) a = 1 and ( 4 ) a 7^ 0 , 1 . 

( 1 ) Consider when a = 0 . After splitting with respect to it can be 
verihed that 021 = 0 and one is left with the following determining 
equations 


((1 - 27 )(Ciw + an))/ - a^ug = 0, 

2(7^11 + (7 - 2)6)/ + (^22 + 6 + C2u)uf = 0, 

7(^12 + Ciu)g = 0 , 

(27011 + (27 - 3)6 + a22)g + (^22 + 6 + C2u)ug' = 0. 


( 24 ) 


From the third equation, notice that if 7 = 0 , G becomes a func¬ 
tion solely of 2; and hence, this case is reducible. Thus, it follows that 
012 = 0 and 6 = 0 . Dividing the second equation by f' and u, and 
differentiating it with respect to m 2 times, one can study the following 
cases: 

/ 


(a) 

(b) 


7^ 0 and 


'x: 

For the case when 


= 0 . 


f \" (2 — 7) 

-' 7^ 0, it follows that on = 6"^- 


\uf'j 7 

Consequently, C2 = 0 and 022 = —6- Substituting this into the 

determining equations ( 1241 ) . no other extensions of the generator 
are obtained apart from the studied subalgebra. 

^ f \" r. ..... .. . f 

= KU + p. 


For the case when 


= 0, it follows that 


uf'J uf 

Furthermore, the general solution of this depends on f 5 . Thus, one 

needs to study whether 7^ 0 or = 0. 
i) For the case when 7^ 0 , the general equation for / (with 

fiX 

a possible shift) is /o — . Substituting this into the deter- 


,uj 

mining equations, one gets 022 = 


(27 - /3 - 4)6 + 27011 
/3 


and 
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( 2 ) 


C2 — 


2/v(7aii + (7 - 2)^1) 


Consequently, the general solution 


/1\/ 3+1 

for is 5^0 ~ • The extension I 5 X^ + 2 ^Xq is obtained along 

\uj 

with the studied subalgebra. 

ii) For the case when = 0 , it follows that /t 7^ 0 . Hence, the 
general equation for / is Substituting this into the de¬ 


termining equations, one obtains C2 = 


2(7011 + (7 - 2 ) 6 ) 


and 


K, 

022 = Consequently, the general solution for g is goe'^^'^. 
The extension k,X^ + 27X4 is obtained aside the studied subal¬ 
gebra. 

For the case when cr = ^, after splitting with respect to y, it follows 
that 021 = 0. Also, since (1 — 4 :j)g + ug' = 0 leads to a reducible case 
it then follows that 6 = 0 - The remaining determining equations are 


2012(1 - 27)/ - 2oi2W5' + (012 - 2C,2U^)uf' = 0, 

4((2 - 7)6 - 7011)/ + (oil - 2 o 22 - 6)w/' = 0, 

(1 - 47)0125^ + (012 - 2 C , 2 U ^) ug ' = 0 , 

((1 - 47)011 + (7 - 47)6 - 2 a22)g + (on - 2022 - ^i)ug' = 0 . 

Dividing the fourth equation by g (as it is nonzero) and differentiating 
it with respect to u, one is left to study the following cases: 

(a) {X\ ^Oand 

a, 

j 7^ 0 . It follows that on = 2022 + 6- H 7 = 0 

/I — 

then 6 = 0 , but if 7 7^ 0 then 022 = 6 - • From the third 

equation, one needs to study the following cases: 

(i) the case where there exists a generator with 012 7^ 0, and 

(ii) the case for which all generators have 012 = 0. 

If there exists a generator with 012 7^ 0 , then g satishes 

the form (1 — 47)5^ + (1 — ( 5 v?)ug' = 0 . Notice that (5 = {) 


Consider 


9 
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( 3 ) 


is reducible. Hence, ( 5 ^ 0 . Without loss of generality, one 
can assume that (5 = 1 . Then the general solution of g is 


9 q 


1 - 


1 


where 7 = 


1 — 47 


7^ 0 (if 7 = 0, the case is re- 


j 2 

ducible). Substituting this into the determining equations, we 

1 \ 7 +( 1 / 2 ) 


^12 


obtain C2 = —. It follows that / = 0 (n) 1 — 


where 




— fo — ‘^90 



Here, the extension Xa + 2X7 is 


l)(l/2) 

obtained besides the studied subalgebra. 

For the case where all generators have ai2 = 0 , it follows that 
(2 = 0 . No other extensions are obtained for this case. 


Consider 


9 


= 0 . The general solution for this is g = gou" 


Substituting this into the determining equations, one obtains that 
ai2 = 0 and C2 = 0 . Consequently, the form of / either satishes 
+ 1 )/ — uf = 0 or not. If it is satisfied, then the general solu¬ 
tion is / = Moreover, 022 = (/t — dq + l)(aii — ^1) + 8^1. 

Here, the extension [k, + 1 )X2 + 2X5 is obtained apart from the 
studied subalgebra. If / does not satisfy {k + 1 )/ — /'u = 0 , then 
no extensions are obtained other than the studied subalgebra. 

For the case when cr = 1 , the determining equations after splitting 
with respect to y are as follows 


(1 - 27)Ci/ + (Cl - C2u)uf' = 0, 

((1 - 27)012 + ((4 - 27)^114 - 2701114))/ - 012145 ^ 

+((oii - 022)^ + ai2 - 02114^)14/' = 0 , 

(1 - 27)Ci 5' + (Cl - C2u)ug' = 0, 

- 02114 / + g{{l - 27)01114 + (1 - 27)012 + (4 - 27)/ i14 - 022^) 

+5''i4((oii - 022)w + ai2 - a2iu‘^) = 0. 

From the first and third equations, one can study the following 2 cases: 

(a) fg' - gf = 0, and 

(b) fg' - gf' 0 . 

Notice that when fg' — gf = 0 , then g = gof. This is a reducible 
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case. Hence, we consider only when fg' — gf ^ 0 . It follows that 
Cl = C2 = 0. From here, one can assume that g = 0(w)/ (as / is 
nonzero), where (f)' 7^ 0 . If it is assumed further that 0 = i/^(w) + l/w, 
then the determining equations are reduced as follows: 


(2(-7aiiw + (2 - 7)Ciw) - (7 + 1^11)012)/ 

+((aii - 022)11 + ai 2 - a 2 iu^)uf' = 0 , 

(oil - 022)11 + ai2 - 02iii^)iC' + oi2iC^ + (oii + 2oi2ii“^ - 022)1/^ = 0. 

( 25 ) 

These equations lead one to study the two cases where: 

(i) there exists at least one generator with 012 7^ 0, and 
ii) where all generators have 012 = 0. 

For the case where there exists at least one generator with 

+ Xu + (3 


012 7^ 0, it follows that iC(ii) = 


where (3 7^ 


u{P — iCoii) 

0 , 1/^0 7^ 0,X, K are constant. Without loss of generality, it is 
assumed further that 13 = 1 . Consequently, we obtain on = 
A012+022 and 021 = —K,ai2. Substituting this into the remaining 
determining equations gives the solution for / which depends 
on the following three cases: 

(A) 4 /v — > 0 , 

(B) 4 k, — < 0 , and 

(C) 4k-X‘^ = 0 . 

For the case where 4 /v — A^ > 0 , it is assumed that 4 k — 
A^ = p 7^ 0 . The solution for / is 


(1 - ^ 

^ {ku‘^ + Aw + 1)"^ 


(2A7- 4 ^) 


arctan 


P 


' X + 2ku'' 
, P / 


where g is constant. 

For the case where 4 w — A^ < 0 , it is assumed that 4k — 
A^ = —p^, p 7^ 0 . The solution for / is 


/o 


(1 — ^ f 2ku 4 - X — p 


P 


{ku^ + Aw + l)"^ \2(3ku + A + p, 
where /io is constant. 
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For the case where 4 av — = 0 , it follows that 

/ 4(7 + /iM) \ 

(1 - 'il)ou)u^^-^ y Xu + 2 ) 

^ + Aw + 1 )^ 

where /i is constant. 

If 7 7^ 0 , then 022 = - -If 7 = 0 , then = 

7 

——. Here, the extension (A7 — fi)X^ — /iXq + 7X7 — 

Zu 

is obtained apart from the studied subalgebra. 

For the case where all generators have ai2 = 0 , the determin¬ 
ing equations are reduced to 


2((2 - 7)^1 - 7 aii)/ + (ail - ^22 - a2iu)uf' = 0 , 

((4 - 27)^1 - 022 + (1 - 27)011)5' - 021/ ( 26 ) 

+(aii - 022 - a2iu)ug' = 0. 

Dividing the hrst equation ([2SD by uf' and differentiating this 
equation with respect to u twice, leads to the study of the 
following sub-cases: 

/ 


(A) 

(B) 


.ufX 

'x: 

If I 


7^ 0, and 


= 0 . 

/ 


^uf 
2-7 

6 - 5 ^22 = 6 


7^ 0, then it follows that if 7 7^ 0 then on = 

2-7 


and 021 = 0 . If 7 = 0 then ^1 = 0 , 


7 7 

022 = Oil and 021 = 0 . For both cases, no extensions are 

obtained apart from the studied subalgebra. 

/ / \" 

= 0, then the general solution for / is /o 


If 


u 


J'u, 


.1 + u, 


where 7^ 0 (else it is reducible). Substituting this into 
the determining equations ( 1261) . one obtains that 021 = 


-7011 + (2 - 7)^1 




and 022 = 


[ k , - 27)011 + (4 - 27)^1 




Substituting this into the remaining determining equation. 
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one finds that g satisfies g'u{l + m) + (1 — n)g + / = 0. The 


general solution of this is g = ^5^0 “ /o 


u 


u 


n—1 


u ~\~ \) VmTI/ 


The extension kX2 + 2 {Xq + X^) is obtained aside from the 
studied subalgebra. 

( 4 ) For the case where a 7^ 0 , 1 , the determining equations are split 

with respect to y. Since /' 7^ 0 , it follows that C2 = 0 and 021 = 0 . 
Notice also that since aug' + {a — 2 'y)g = 0 leads to a degenerate case, 
then Cl = 0 and ai2 = 0 . The remaining determining equations become 


{aau + (a - l)Ci - ^22)^/' + (-27011 + (4 - 27)^1)/ = 0 , 

{aan + {a - l)Ci - a22)ug' + ((a - 27)011 + (a - 27 + 3 )Ci - 022)5^ = 0 . 

. . 

Dividing the first equation by / (as it is nonzero) and differentiating 
with respect to u, it can be observed that there is need to study the 
following cases: 


(a) 

(b) 


'uf 

f 

'uf 

f 


7^ 0 and 


= 0 . 


For the case with 


V'' 

/ 


7^ 0, it follows that 022 = <aoii + {a — 


l)Ci. Substituting this into the remaining determining equations 
_ 2-7 

(IzZp . we find that if 7 7^ 0 then on = -Ci? if 7 = 0 then 

7 

Cl = 0 . Substituting all these, no other extensions of the generator 
is found other than the studied subalgebra. 

'uf' 


For the case where 


/ 


= 0, the general solution for / is 


where k ^ 0 . Substituting this function into the determining equa- 

,. ,, , (/ta - 27)011 + (/ta -- 27 + 4 )Ci 

tions, it follows that 022 = - 


K 


This leads us to study the two cases, that is, if g'u — g{K — 1 ) = 0 
or g'u — g{K, — 1 ) 7^ 0 . If g'u — g{K — 1 ) = 0 , the general solution 
for g is gou'^~^. Another extension of the generator apart from the 
studied subalgebra is found, that is KX2 + 2 Xq. For the case where 
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g'u — g{K — 1) 7^ 0, if 7 7^ 0 then an 


^7 

7 


^1, and if 7 = 0 


then = 0 . Again, these lead to a generator with only the studied 
subalgebra as its extension. 


The complete representative classes for the autonomous system with all 
admitted generators having = 0 is listed in Tables [ 2 ] and El 


4 Conclusion 

A complete group classihcation of the systems of two autonomous non¬ 
linear second-order ordinary differential equations of the form y” = F(y) 
excluding the systems which are equivalent to linear systems and the de¬ 
generate case were presented using both the direct and algebraic approach. 
The important thing in this study is that the analysis of the determining 
equations where split into two cases: 1) the case where at least one admit¬ 
ted generator has 7^ 0 and 2) the case where all admitted generators 
have = 0 . The hrst was analyzed through the direct approach while the 
latter was analyzed using one-dimensional optimal system of subalgebras 
followed by the direct approach. For the direct approach, all possible Lie al¬ 
gebras were found with the aid of the equivalence transformations applied 
on the determining equations. As for the algebraic approach, the study 
was reduced to the analysis of relations between constants of the generator 
with its corresponding basis operators. The obtained classihcation is sum¬ 
marized on Tables [U [ 2 ] and El It is highly likely that the same methods 
shown in this paper are applicable to the group classihcation of systems of 
two nonlinear second-order ordinary diherential equations, which will be 
next goal for further studies. In addition, it is also believed that this can 
be extended to systems in more general cases. 
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Table 1 

Group classification of systems admitting at least one generator with 7^ 0 . 


F 

G 

K 

Extension of Kernel 


Kz + 

Z'^ \ y / 

0 

T3, y4 



-1 

Yj, Ys, n 



1 

>"9, no, n 

i^y + {foy - fiz)T{y, z) 

Kz + {foz + /iy)T(y, z) 

0 

n, n, n 



-1 

n, n, n 



1 

n, no, n 

Ky + e~- z~^{foy + fiz) 

Kz + foz~^e^ 

0 

n, n, n 



-1 

n, n, n 



1 

n, no, n 
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Table 2 

Group classification of systems admitting all generator with = 0 . Here we have 6 i{u,v) = (cos(n)/(n) + sm{u)g{v)), 


Oi i 

6'2(n,n) = sin(n)/(n) - cos{u)g{v), xi(a) = o , i ^2(a) = o , . • 
_ + 1 _+ 1 


F 

G 

Relations 

Extension of Kernel 



u = ^ —1<q;<1 

2 — — 

7X2 + X5 + aXg 



u = 2/e“^ 

7X2 + X4 + X5 

v) 

—n) 

y = ncos(n), z = nsin(n) 

7X2 — X7 + Xg 

e(“-2T')“0i(n,n) 

g(“-27)«02(u, n) 

y = t!e"“ cos(n) + xi (a) 

2; = ne"“ sin(u) — X2(o:), a > 0 

7X2 - X3 + a(X5 + Xe) - X7 + Xg 

e(“-2T')“0i(n,n) 

e(“-27)«02(n, n) 

y = t!e"“ cos(tt) — xi (ck) 

2; = ne"“ sin(n) + X2(o:), a > 0 

7X2 + X3 + a(X5 + Xe) - X7 + Xg 

e(“-2T')“0i(n,n) 

g(a- 27 )n 02 (u,n) 

y = ne"“ cos(n) 

z = ne"“ sin(tt), a > 0 

7X2 + a(X 5 +X 6 )-X 7 + X8 

{g{v)u + 


y = nn, z = n 

7X2 + X7 

+ /(n))e(- 272 ) 


u = — 2 y 

7X2 + X4 + X7 

{{y/z)g{u) + /(u))e((^"27)(yA)) 


u = ze~y/^ 

7X2 + X5 + Xq + X7 


g{z)e-^'^y 


7X2 + X3 



Table 3: Group classification of systems admitting all generator with = 0. Here, /o, go, (po, 4>i, ca, (3, k, go-, ^ and 7 are constant. 


F 

fozPy^+^ 


Subalgebra 1 . 7X2 + 

G Relations Additional Extension of Kernel 

goz^+^'-y'^ 7 = -27 ^ 0, /3 5^ 0 - 7X5 

7 = -27 5 ^ 0, K ^ 0 kXs - 7 X 4 


F 

4>{u){y - z'^y 


/02-(«+l)y^+l 


Subalgebra 1 . 7A2 + 




G 

9o{y - 

goz-^y^ 


7 = 


Relations 
1 — 47 


7^0 


4 ' = fo{y - z^y^^ + 2goz 
, fc + 1 — 47 , 

/^ + 1 7^ 0, 7 = - 7^ 0 


Additional Extension of Kernel 

X 4 + 2 X 7 

{k + 1)A'2 + 2Xq 


to 

00 


z-ay 

° (22 + \yz + Ky^)-t 




-27 


-fo 


Subalgebra 1 . 7A2 + X 5 + Xq 

G Relations 

{Ky + (A + 0 ) 2 ) 


( 2 ^ + Xyz + Ky'^)'< 


Yi{y,z) 

2 A 7 — 4/i A 2 + 2/^y 

-arctan- 


2 = 1,2,3, o 7 ^ 0, 


Here, ipi {y, z) = e P 


' 42 {y,z) = 


2 /^ —A 7 

2Ky + {X + p)z'^ p 
2ny + (A - p)zj 


with 4 k — A2 = p2 p ^ g. 
with 4 k — A2 = —p2^ p ^ 0; and 


4(mP + 7^) 

■ipsiy, 2 ) = e Ap + 22 with 4 k — A^ = 0 




7 5^ 0, K ^ 0 


Additional Extension of Kernel 
(A 7 — p,)Xs — p,Xe + 7 X 7 — K'yXs 


kX2 +2(X6+X8) 


F 

/o2-«p^+i 


G 

goz^-'^y^ 


Subalgebra 1 . 7A2 + A5 + aAg, o; 7^ 0 

Relations 

7 = cuK — 27 , a 7 ^ 0 , 1 / 2 , 1 , 7 ^ 0 


Additional Extension of Kernel 
1^X2 + 2 Xq 

Table 3 - continued on next page 



Table 3 — continued from previous page 


F 


Subalgebra 2 . 7X2 + X4 + X5 


G 


Relations 
a. — K 


9 oy^ 


7 = 


- , KO. 7^ 0 


Additional Extension of Kernel 
0X5 + /^A’4 


F 

(/o cos(ti) + g<3 sin(n))ij'' 


(/o sin(n) 


Subalgebra 3 . —JA7 + As 
G Relations 

- go cos{u))v'-^ fo go^ 0, 

u = arctan(2:/y), 


Additional Extension of Kernel 

—-—A2 + A5 + Aq 

__ 


to 

CO 


e^^(/o cos(ri) + ^0 sin(ri))'i;' 


-7K —3 


Subalgebra 3 . 7A2 — A7 + As, 7 7^ 0 
G Relations 

e^“(/o sin(M) - 30 cos(m))i ;-'^«-3 /o ^ 0, 30 ^ 0, 

7 = —27 7^ 0, = arctan(2/y), + 2;^ 


Additional Extension of Kernel 

2A2 + As + Ae + ^(As - A7) 


Subalgebra 4 . 7A2 — A3 + q:(A 5 + Ag) — A7 + As, a > 0 
F G Relations 

cos(u) + go sin(w))i;'^ sin(w) — go cos(u))v*^ /o / 0, ^'o 7^ 0^ 


Additional Extension of Kernel 


u = arctan 


/ z + X2(q) \ 

V3-Xl(o)/’ 


= e-2“'“((3 - XI(o))^ + (2 + X2(«))^), 
Xl(«) = 9°", . X2(a) = 


1 — . 


1 ' 


+ 1 


- A 2 + As + As — Xl^3 + X2^4 


Subalgebra 4 . 7A2 + A3 + a(As + Ag) — A7 + As, a > 0 

Table 3 - continued on next page 



Table 3 — continued from previous page 


F G 

g(o'-27)-u^y^ cos('u) + QQ sin(w))i;'^ sin(«) — cos('u))tj” 


u = arctan 


Relations 

/o 5 ^ 0, go 0, 

^ - X2(a)' 


Additional Extension of Kernel 


/ z-X 2 ia) \ 

Vg + Xl(a)/’ 


-.y + xiio‘)^ 

+ XI («))^ + (2 - X2(«))^), 
a , , 1 


Xl(“) = 




1 


X2(a) = 


1 — K . 


+ 1 


-X 2 + Xs + Xq + Xl^3 — X2^4 


F 

^{a-2'y)u cos(u) + QQ sin(w))i;'^ 


Subalgebra 4 . 7A2 + 0(^5 + Ae) — X'j + Ag, o > 0 
G Relations 

e(““^^)“(/o sin(M) - go cos(«))i;'^ /o ^ 0 , go 0 , 

u = arctan(2:/y), + 2^) 


Additional Extension of Kernel 

—-—A2 + A5 + Ag 

__ 


00 

O 


F 

gQZ^~^e~‘^^^{y + k^z) 


Subalgebra 5 . 7A2 + A7 


G 

goz^e-y^^ 


Relations 
7 = 27 


Additional Extension of Kernel 

As + 7^6 + (/3 - 1)^7 


F 

(go2 + /o)e(^“-2^") 


G 

goeW^-^T-^) 


Subalgebra 5 . 7A2 + A4 + A7 

Relations 

u = -2y, 13^0 


Additional Extension of Kernel 
/3X2 + X3 


F 

(ff 02 + foils + uy/^)(/3 + u}'‘ 


G 

goiP + uT 


Subalgebra 5 . A4 + A7 

Relations 

u = z^ — 2y, K ^ 0 


Additional Extension of Kernel 

(1 - 2a.) A2 + 2 i-i 3 Xs + 2A5 + Ae) 


Subalgebra 6. 7A2 + A5 + Xq + A7 

Table 3 - continued on next page 



Table 3 — continued from previous page 


F 

(goy + /oz)^”' 


G 

goz'^e-y(v/^'> 

Relations 

7 = 27 + fc — It^O 

Additional Extension of Kernel 

(k- 1 )A 2 - 2 (A 5 +A 6 ) 


F 

){kz + 70 i) 

-^ymz) 

G 

gozPe’^’^-'^y 

g^^Pz + I^Z^-jy 

Subalgebra 7 . 7X2 + Xz, 7 0 

Relations 

7 = 27, fo = 90/7 

7 = 27, 

(j> = 4 > 0 Z + <j>i, (j >0 0, K = ( 70 o )/2 

Additional Extension of Kernel 

(13 - 1 )X 3 + nXr + 7^6 

/3X3 + 2 KXr + 7X4 



Appendix 1. For = 0, the determining equations have the form 


I Cl = FyCl + FzC 2 + 9 l, 

I C2 = ^J/Cl + GzC2 + ^2, 

where qi and q 2 are functions of y and Differentiating them with respect 
to X, one has 

I cr = FyCi + 

\c^^ = GyC[ + GzC^, 

Differentiating the latter equations with respect to y and 2 ; 


FyyCl + FyzC2 ~ 0) 
FyzCl + FzzC 2 — 0) 

GyyCl + GyzC2 — O5 
GyzCl + GzzC 2 — O5 


F‘^ F F 

/■' p _ yz _ n ^ — 0 — n ^y^ 

Si) ^yy '^yy '^yz ^ '^yz '^zz ^ 

zz zz zz zz 


Case 1. Let Fzz 7 ^ 0, then 

CH 

Thus, 


^yz >./ 


= 0 . 


F 


yz 


F; 


= k, 


zz 


and 


or 


Fyy kFyz 0 , Gyy kGyz 0 , Gyz kGzz 0 


(Fy — kFz)y — 0 , (Fy — kFz)z — 0 , {Gy — kGz)y — 0 , {Gy — kGz)z — 0 , 


One has 


Fy - kFz = ki, Gy - kGz = /c2, 

dy dz dF 

T 


F = ^{z + ky) + kiy, 


—k ki 

dy dz dG ^ ^ , , x , 

— = — = G = ^{z + ky) F k 2 y. 
1 —k k2 
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Changing the variables 


the original system 


becomes 


y = y, z = z + ky, 


y' = F{y,z), z'' = G{y,z), 


y" = + kiy, z = {k^{z) + "^{z)) + {kki + k 2 )y 

Thus one needs to study the equations 

y” = kiy + F{z), z!' = k 2 y + G{z), k 2 F" ^ 0. 

Cf = ^iCi + F'C;2 + Qij 
C'i = ^2Cl + G'C2 + ^2, 

Notice that because of F" 7 ^ 0, then C 2 = 0 from the second equation 

0 = /C2C1 Cl = 0 . 


Case 2. Let Fzz = 0, then by symmetry Gyy = 0. Hence 

FyyC'l + FyzCl — 0) FyzCl — 0) GyzC2 — O 5 GyzCl ^ G zzCl 

If 

Fyz^^=> C,1 = 0 , C2 = 0 - 


= 0 , 


Hence, 

and 

Thus 


Fyz — 0 , Fzz — 0 , Gyy — 0 , Gyz — 0 


FyyCl — 0 , GzzC2 — 0 - 


y" = kiz + F{y), z" = k 2 y + G{z), kMF"^ + G"^) ^ 0. 

Cl = ^^Ci + ^iC2 + 91, 

C2 = ^sCl + G'C2 + ^2, 
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Let F” ^ 0, then Ci = 0 because 

C = F'Ci + A:iC 2 =» /ciCH 0 ^ CH 0 

Let G" ^ 0, then C 2 = 0 because 

cr = hc'i + G'c^ ^ /C2c; = o ^ c; = o 
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